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Entanglement in multipartite quantum systems, due to 
its importance in quantum information theory, became a 
topic of numerous investigations (see [l| for a recent com- 
prehensive review of various aspects of characterization 
and applications of multipartite entangled states). Until 
now, however, the efforts did not bring final solutions to 
many fundamental problems. One of the question is the 
classification of states which are interconvertible by lo- 
cal unitary transformations, i.e. operations on the whole 
systems composed from unitary actions (purely quantum 
evolutions) each of which is restricted to a single subsys- 
tem. To appreciate the importance of such an experi- 
mental setting let us remaind that it is a basis for such 
spectacular applications of quantum information tech- 
nologies like teieportation or dense codding where the 
fundamental parts of experiments consist of manipula- 
tions restricted to parts of the whole systems in distant 
laboratories. 

Recently, the problem of checking local unitary (LU) 
equivalence of qubit systems was solved in 0,0]. In this 
letter we present a geometric approach to the problem 
which supersedes previous ones by its applicability to 
systems of arbitrary dimensions treating on equal foot- 
ing qubits and qudits. The same geometric formulations 
allows also to answer in arbitrary dimensions the ques- 
tion of distinguishing multipartite states by local mea- 
surements (i.e. measurements constrained to single sub- 
systems) . By taking advantage of relations between both 
concepts we are able to give an universal algorithm for 
checking LU-equivalence of pure states in arbitrary di- 
mensions. Despite a rather abstract geometrical picture 
underlying our approach the final algorithms for check- 
ing equivalence and distinguishability are easy in appli- 
cations. 

Consider a quantum system consisting of M identical 
subsystems with the Hilbert space T-L = ^if^^i where 
each Hi is .ZV-dimensional dimensions of the subsystems 
can be done along the same lines, albeit with . We say 
that two pure states \vi) £ W and \v 2 ) <E H are locally 
unitary equivalent (LU-equivalent) if and only if there 
exist unitary operators Ui € U(iV) such that 

= U 1 ®...®U M \v 2 ). (1) 
We call states |i>i) and \v 2 ) indistinguishable by local 



measurements if and only if 

= (v 2 \Av 2 ) VA (2) 

where A is a local hermitian operator, i.e, it is real com- 
binations of A\ g§ I ® . . . ® /, ... , 1 65 . . . ® I (8 Am and 
Ai are N x N hermitian matrices. Each operator corre- 
spond to a measurement performed on a single subsys- 
tem. Physically, we may imagine that subsystems are 
located in different laboratories, each local measurement 
is performed on a separate copy of the state in question 
and by pooling results obtained from such measurement 
we want to distinguish different states. 

Pure states are, in fact, points in the projective space 
P(%) rather than vectors in H - two vectors \v), \w) £ H 
represent the same state if and only if \v) = z\w), where 
z € C. Identification of all vectors differing by a mul- 
tiplicative complex constant gives a point in P(H). We 
will denote by 7r(|w)) the point in the projective space 
corresponding to \v) 6 H, i.e. it denotes the canonical 
projection from H to P(H). The action of the unitary 
group on % as in (JlJ translates via 7r to an action on 
f(Ti). We can thus reformulate the definition of LU- 
equivalent states in the language of the projective space. 
Two pure states \v\) and \v 2 ) are locally unitary equiv- 
alent (LU-equivalent) if and only if there exist special 
unitary operators Ui <E SU(7V) such that 

Tr(\v 1 ))=ir(U 1 ®...®U M \v 2 )). (3) 

In more technical terms LU-equivalent states lie on the 
same orbit of G = SU(iV) <g> . . . ® SU(AT) action on ¥(H) 
and checking the LU-equivalence can be reduced to in- 
vestigations of the structure of such an orbits. We re- 
stricted the group to SU(N) since vectors differing by a 
phase factor are indistinguishable on the projective space 
level. Furthermore, two states represented by arbitrary 
(not necessary norm one) vectors \vi) and \v 2 ) are indis- 
tinguishable by local measurements if and only if they 

fulfill (GO) with k) substituted by } Vi ' . 

The study of orbits of the group G in the projective 
space F(H) is facilitated by the existence of some addi- 
tional structures on P('H) and G which we now briefly 
describe. 

The projective space P("H) is a symplectic manifold. 
It means that there exists a differential two-form w on 



2 



¥(H) which is closed (duj — 0) and nondegenerate. To 
define lo we have to determine its action on a pair of 
arbitrary vectors A x , B x tangent to P(H) in an arbitrary 
point x. Now, each tangent vector A x is a vector tangent 
to the curve t i-> 7r(exp(£A)|v)) at t = 0, where A is 
some element of the Lie algebra of the full unitary group 
U('H) [in our case by a choice of a a basis in Hi it can 
be identified with U(iV M )], and \v) is a vector in H such 
that 7r(|i;)) = x. We have then 



u(A x ,B x ) 



-Im 



(u|t>) 



2 (d^ 



(4) 



It is easy to check that the above definition is correct, 
i.e. it does not depend on the choice of \v) and more- 
over u) is invariant with respect to the action of [/(%) on 
P('H), i.e. the action of U(H) is symplectic. This fact is 
a basis for another canonical construction in symplectic 
geometry, namely the moment map Let us denote 
bu u('H) the Lie algebra of the unitary group U('H), i.e. 
the algebra of all anti-Hcrmitian operators acting on % 
and by u*(H) the space dual to u(H), i.e. the space of 
all linear operators on u('H). The latter can be identified 
with the space of all Hcrmitian operators (observablcs) 
on % and the action of 7 6 u*(H) on I £ u(W) is given 
by (Y,X) = Tr(XY), where in order to calculate the 
trace we may choose an arbitrary basis in T-L and treat 
all operators as appropriate matrices. 

Instead of giving the formal definition of the moment 
map fi in the most general setting (see we write ex- 
plicitly an expression valid in the considered case where 
/j : P(H) —> u*(T-L) can be shown to read as 



(fi(x),A) 



i (v\Av) 



n(\v)). 



(5) 



Besides of its fundamental mathematical origin the 
idea of moment map has a natural physical interpreta- 
tion. Namely to any state \v) we assign a linear functional 
/j,(tt(\v))) G u*("H) which encodes information about all 
expectation values of all observables in the state \v). In 
the following we will be interested in the restrictions 
of u> and fi to an orbit G.x of the group G through 
a point x G P(H). They can be calculated as in (0| 
and ([5]) but now A and B are restricted to elements of 
q = su(7V)©. . .©su(iV), i.e., the Lie algebra of G. In case 
of the moment map it means that to every state \v) we 
assign an element of g* which encodes expectation values 
of Hcrmitian operators of the following type 



Yx®I<. 



I 
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(0) 



It means we extract from state \v) the information about 
local measurements as we only control operators of ([6]) 
type. 

Observe that we can identify two natural actions of 
the group U{H). The first one is the action on P('H) 
induced form the ordinary action on % itself, whereas 



the second is the coadjoint action on u*(H) defined as 
Y i — ^ UYUK Both are intcrwined by the moment map. 
It can be directly checked using ([5]) that ji is equivariant, 



(n(*(U\v))),A) = (Un(i r (\v))U*,A). 



(7) 



The moment map projects thus the orbit through x = 
7r(|i>)) on the orbit of the coadjoint action through 
iu(7r(|i>})). As already mentioned we may restrict the rea- 
soning to the actions of G which is a subgroup of U(H). 
Generally the mapping by the moment map will not be 
a diffcomorphism between the orbit G.x and the coad- 
joint orbit through fi(jr(\v))). The latter operator can be 
stabilized by some subgroup Stab([i(ir(\v)))) of G, i.e., 
Ufj,(ir(\v}))W = ^(7r(|t;})) for some nontrivial U. Like- 
wise, we may have 7r(y|w)) = 7r(|w)) for some nontrivial 
unitaries V forming a subgroup Stab(ir(\v))) of G. In 
general the stabilizers Stab(ir(\v))) and Stab(fj,(ir(\v}))) 
are not equal, so the mapping is not diffcomorphic. The 
other facet of this fact is that the form uj restricted to 
G.x is no longer nondegenarate; G.x is not a symplectic 
space. The difference between dimensions of both stabi- 
lizer subgroups conveys an important information useful 
in deciding the LU-equivalence as will be shown in the 
second example further in the text. 

The fact that /i maps orbits of G in H on orbits of 
the coadjoint action is crucial for identifying states on 
the same orbit. The relevant observation is that each 
coadjoint orbit intersects the subspace in g* which is dual 
to the maximal commutative subalgebra of g In our 
case this subspace consists of operators of the form (|6]) 
with Yi diagonal. Let us explain how this fact is reflected 
in the present setting. To this end we consider a state 



N 



H, -,»M=1 
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^From the coefficients C^...^ we can build M following 
N x N Hcrmitian positive scmidefinite matrices 



ikjk 



G G 

^l\...% k ...l M ■ il-.-Jk-.-lM ' 



(9) 



where by the overbar we denoted the complex conju- 
gation and the summation over repeating indices is as- 
sumed. The orbit of G through x = k(\v)) is now mapped 
by fi on some coadjoint orbit and the fact that the latter 
contains a point ([6]) with Yi diagonal means, after trans- 
lating back to the orbit G.x, that it contains a point 



i M ) for 



, I<M=1 J*i>®.. 
which the corresponding matrices C' k are diagonal, C' k = 

diag(j^ fc , . . . ,p 2 Nk ). The numbers pf k have a straight- 
forward interpretation as probabilities to find the fc-th 
subsystem in the state \l) while neglecting other M — 1 
subsystems (see [H for detailed calculations). Let us re- 
mark that in the case of two subsystems, M = 2, the 
whole procedure leads to the familiar Schmidt decompo- 
sition which allows for a straightforward identification of 



3 



LU-cquivalcnt states for bipartite systems. There is no 
useful generalization of the Schmidt decomposition for 
more than two subsystems, nevertheless, the following 
reasoning shows how to overcome this obstacle. 

To pass from \v) to \v') one has to act on \v) by 
U-f (8 . . . ® U M ( T means the transposition), where Uk £ 
SXJ(N) are such that UlC k U k = diag(p 2 lk , . . . ,p 2 Nk ). It 
is always possible to find such unitaries since, as already 
mentioned, Ck are Hcrmitian. Moreover one can always 
choose Uk in a such a way that p\ k > . . . > p 2 Nk and 
in the following we assume it was done. In fact all in- 
tersections of the coadjoint orbit with © with diagonal 
Yi differ only by ordering of the diagonal elements. We 
avoid thus this ambiguity by fixing the order. 

The image of tt{\v')) under the action of the moment 
map n can be calculated @ . It has the form © with 



Yfc = adiag(- 



N 



1 

N 



' PNk)j 



(10) 



where a is some appropriate const. Going back to the 
above mentioned physical interpretation of fi we see that 
in fact moment map encodes information about the state 
\v') contained in local measurements as it is determined 
only by the local probabilities p kl . 

The matrices C' k and, consequently Yk have, in general, 
degenerate spectra, i.e, several pik can repeat in (|10[) . Let 
us denote by vu the number of different eigenvalues of C' k 
and by m,k,n the multiplicity of the n-th eigenvalue. It is 
now easy to show [f| that the stabilizer Stab{ji{-K{\v')))) 
consists of U\ ® . . . (8 Um € G where 



Uk.o 



Uk,l 



\ 



Uk„ 



\ 



J 



(11) 



with Uk, n S U(mfc : „). The equivariance of [i implies 
that Stab(ir(\v'))) C Stab((i(ir(\i/)))). Consequently, for 
any \v) the corresponding state \v'} is given up to the 
action of Stab(fj,(ir(\v'}))), which is explicitly known. 

We thus arrive at the following algorithm to check if 
two states |i>i) and |t>2) are LU-equivalent. 

• For both states compute the matrices Ck © and 
check whether they have pairwise the same ordered 
spectra. If this is not the case the states are not 
LU-equivalent. 

• If the first condition is fulfilled use the unitary ma- 
trices Uk diagonalizing Ck for each state to find 
\v[) and \v' 2 ). If they are equal then and 
|w2) are LU-equivalent, otherwise they are LU- 
equivalent if and only if there exist U\ <8 . . . <8 Um £ 
Stabiniirdv't)))) such that = Ui <g). . .®U M \vi). 

In case of generic states Stab(fi(n(\v' i }))) consists of di- 
agonal matrices diag(e Ic ^ 1 , . . . , e l ^ N ) and above method 



can be easily applied. For other states our method also 
simplifies computations since we no longer need to use 
the full U(N) (8 . . . <8 U(N) group but can restrict our 
attention to an appropriate subgroup. The worst case is 
when the state \v) is such that all Ck are ^=J- Then the 
action of U(N) (8 . . . (8 U (N) on such state gives a state 
for which again all Ck are normalized identity matrices. 

Our approach gives also a beautiful mathematical char- 
acterization of states which are indistinguishable by local 
measurements. As we noticed before for any state \v) the 
moment map fj, assigns an element of q* which encodes 
all expectation values of local measurements in state \v). 
In other words all states that have the same expectation 
values for local measurements are sent by the moment 
map to the same point of q* - they constitute a fiber over 
this point. To examine if states \v[) and \v' 2 ) obtained 
in the second step of the above algorithm are indistin- 
guishable by local measurements it is enough to check 
if ordered spectra of corresponding matrices Ck are the 
same. An intimate connection between LU-equivalence 
and indistinguishability by local measurements appears 
to be helpful in deciding the former in various situations, 
also in the 'worst case' mentioned above, as we show in 
the second example below. 

A detailed analysis showing how to apply our method 
to different kind of states will be published elsewhere 0] . 
Here we give only two simple examples. The first one 
shows how the proposed algorithm works in a generic 
case and involves three qutrits. The Hilbert space is H = 
C 3 (8 C 3 (8 C 3 and has the real dimension dimffi('H) = 
54. The dimension of the projective space F(H) is thus 
dimij(P('H)) = 52. We are interested in the orbits of 
G = SU{3) (8 SU(3) (8 SU{3) - action on ¥{H). The Lie 
algebra q of G is spanned by A (8 / (8 I, I <8 A <E> I, and 
I ® I ® A with A e su(3) and dim R (g) = 24. Let us 
consider two states 

|*) = a\W) + /3|222) = a(|001) + |010) + |100>) + /3|222), 
|$) = \/2a|000) + a|lll) + /3|222), 

where 3|a| 2 + |/3| 2 = 1 and \/3\ < \a\. Direct computations 
show that for both \ty) and |$) we get: 



C\ — C'2 — C3 



2\a\ 2 
\a\ 2 
\B\ 2 



(12) 



Since the corresponding matrices for both states are the 
same, |4>) and I'J) are indistinguishable by local mea- 
surements. Remember that if for at least one pair of the 
matrices Ci calculated for |<E>) and |^) the results differ 
the algorithm ends deciding that the states are not LU- 
equivalent. Since this is not the case we have to invoke 
the second step. The matrices (|T2|) are already diagonal 
so we do not need to change the form of |^) and |<I>). Let 
us thus look at the states that are LU - equivalent with 
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|$) and for which matrices C±, G2, G3 are given by (IT21) . 
Such states are generated by action of Stab(/i(|$})) on 
|<E>). But Stab(//(|$))) consists of matrices U 1 <E)U 2 <E> U 3 
where: 










a i<t>k2 






-i(<Pkl+<l>k2) 



(13) 



The action of ([13)) on |<f>) gives 

|$') = V2ae l ' pl |000) + e l02 a|lll) + e^ 3 /3|222) (14) 

where 0j are appropriate combinations of 4>ki ■ It is easy 
to notice that \ i S>) is not of the (fT4]) form and in the effect 
states |$) and |$) are not LC/ - equivalent. 

Fibers of the moment map consisting, as mentioned 
above, of locally indistinguishable states are given as 
common level sets {Ja(x) = ca} of functions /_a(7t(u)) = 
(/j,(\v)),A) [see ©], where A £ g and c,4 are some con- 
stants. An interesting question which one should ask is 
what is the relationship between the tangent space to the 
fiber of the moment map at given point x = 7r(|u)) and 
the tangent space to the orbit of G-action at the same 
point. These two tangent spaces generate infinitesimal 
movements in P('H) in the direction of states which are in- 
distinguishable by local measurements from state x and, 
respectively, L [/-equivalent with x. The tangent space to 
the fiber at a given point x = tt{\v)) is of course contained 
in the intersection of kernels of all 1-forms d x fA, where 
A is an clement of g. It is very important to emphasize 
that it might be just a subspace of this intersection Q. 

To continue our discussion we will need the following 
definition. 

Definition 1 Let (V, uS) be symplectic vector space and 
W C V its subspace. Then uj- orthogonal complement of 
W is given by 



W J 



{v £ V : oj(w, D) = 0,Viue W}. 



(15) 



In [6[ we prove that the tangent space to the orbit of fr- 
action at x is spanned by vectors which are w-orthogonal 
to the intersection of kernels of all 1-forms d x fA for A £ g. 
Knowing just these kernels we can thus find the dimen- 
sion of an orbit of the G-action. The other facet of the 
above reasoning is that when the tangent space to a fiber 
of the moment map is contained in the tangent space to 
the orbit G.x then the local unitary equivalence is the 
same as the indistinguishability by local measurements, 
i.e., in this case we can restrict our method just to the 
comparison of ordered spectra of Gfc matrices. In [H[ we 



prove that the dimension of the part of the moment map 
which is contained in the orbit of the G-action is equal to 
difference between dimension of stabilizers of 7r(|i;)) and 
fi(ir(\v))) mentioned above. 

The next example shows a practical application of the 
above observation. Let us consider a three-qubit state 



|G#Z) = -^(|000> + |111)) 



(16) 



Here dim R (P(H)) = 14. One calculates d = C 2 = G 3 = 
^12x2- Hence Stab(fi(\GH Z)) is the whole local unitary 
group, i.e, Stab(fx(\GH Z)) = SU{2) <g> SU(2) ® SU(2), 
and its action on \GHZ) is the whole orbit O^Qffzy 
Thus the part of the fiber of the moment map contained 
in the orbit @\Qffz) ^ s exa ctly the orbit @\(jffz)' 
intersection, K, of kernels of dfA®i® 1, dfi®A®i, rf/j®J®A, 
where A £ su(3) is seven dimensional. The dimen- 
sion of the cj-orthogonal complement of K is the di- 
mension of the G-orbit through \GHZ) and is equal 
dim(0|£<^ ) = 7. The maximal possible dimension 
of the fiber is also 7 since it is the dimension of K. This 
means that the whole fiber of the moment map is con- 
tained in 0\QHZ)~ ^ e m ^ er thus that two three-qubit 
states for which G\ = C 2 = G3 = ^1-2x2 are automati- 
cally L [/-equivalent, although the situation corresponds 
the 'worst case' when deciding the LU-equivalence should 
be the most difficult. 

The presented examples do not exhaust the variety 
of possible situations met when checking LU-cquivalence 
0. They correspond to, in a sense, two extremal cases, 
and show clearly that dimensionality of subsystems does 
not play a crucial role in applications, although, admit- 
tedly can make calculations more cumbersome. 
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